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ABSTRACT 



We reconsider the ellipsoidal-collapse model and extend it in two ways: We modify the treatment of the external gravitational shear 
field, introducing a hybrid model in between hnear and non-hnear evolution, and we introduce a virialisation criterion derived from 
the tensor virial theorem to replace the ad-hoc criterion employed so far We compute the collapse parameters and Ay and find 
that they increase with ellipticity e and decrease with prolaticity p. We marginahse them over the appropriate distribution of e and p 
and show the marginalised results as functions of halo mass and viriahsation redshift. While the hybrid model for the external shear 
gives results very similar to those obtained from the non-hnear model, elhpsoidal collapse changes the collapse parameters typically 
by (20 . . . 50)%, in a way increasing with decreasing halo mass and decreasing viriahsation redshift. We quahtatively confirm the 
dependence on mass and viriahsation redshift of a fitting formula for 6^, but find noticeable quantitative differences in particular at 
low mass and high redshift. The derived mass function is in good agreement with mass functions recently proposed in the hterature. 
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1. Introduction 



The 
199: 

Bartelmann et al 



pherical-collap se mo 
Engineer et al.' '20001; 



mode l (e.g. 



2006; 



Mota & van de Brucki 
Schafer & Kovama 2008) is a funda 



Wang & Steinhardf 
l2004t 



mental ingredient in the theory of cosmic structure formation. 
Following the collapse of a slightly overdense, homogeneous 
sphere, it allows the derivation of two essential parameters; 
the overdensity Ay of a virialised halo compared to the mean 
or the critical cosmic density, and the critical linear density 
contrast 6c- The former is important because it allows relating 
sizes to masses of virialised structures, and the latter because 
it establishes a link between linear structure formation and the 
population statistics of collapsed haloes. Despite fundamental 
doubts as to the validity of such a simplified model for accurate 
cosmological predictions, the parameters derived from the 
spherical-collapse model or variants thereof allow surprisingly 
far-reaching predictions such as the halo mass function or 
the correlation properties of haloes, which are confirmed by 
numerical simulations. 

The statistics of a Gaussian random field implies that spher- 
ical collapse should not occur. In fact, the probability distribu- 
l ion for the eigenvalu es At of the Zel'dovich deformation tensor 
dPoroshkevic hl 1 970l) shows that spherical collapse has a vanish- 
ing probability. 



p(Ai,A2,A3) = 



15^ 



SttVSo 



exp - 



36^ 



15/ 
2^ 



(1) 



X (Al - A2){A2 - AiXAi - Ai) 



because it requires Ai - A2 - /I3. In this equation, 6 = Ai + 
A2 + A3, 1 = A1A2 + A2A3 + A1A3, and cr^ denotes the variance of 
the matter power spectrum smoothed on a scale that corresponds 
to a halo mass M. Modifications of the original Press-Sche chter 
mass function dPress & Schechteil[T974t iBond et al.lll99lh mo- 



tivated by ellipsoidal collapse dSheth & TormenlF2002h substan- 
tially improve the agreement of analytic predictions on the halo 
populatio n with numerical sim ulations (e.g. the Millennium sim- 
ulation bv lSpringel et al.ll2005b . 

EUigsoidal collapse was a na lysed many times before 
(see iSi rtelmann et aj Il993t lEisenstein & Loebl 119951: 
Bo nd & MversI 119961 for examples). Several authors have 
worked with the model by iBond & MversI ( 11996) . generalising 
it for different c osmologies and introducing the scale factor a as 
a time variable ('Monaco"1 99l [19971 [l99l ISheth et al.ll200ll; 
Sheth & Tormen 2002; Oht a et al.ll2()04r We are reconsidering 
it here for two reasons. 

First, we want to analyse how different assumptions on the 
treatment of the environment of a halo impact on the parameters 
Ay and 5c- Previous assumptions have been that the principal- 
axis system of the homogeneous ellipsoid is either identical with 
that of the external gravitational shear field, or that the two eigen- 
systems do not coincide, introducing rotation and a deviation 
from the homogeneous mass profile. We introduce another as- 
sumption here, letting the eigensystems of the collapsing ellip- 
soid and its surrounding shear field follow each other until turn- 
around of the major principal axis and then decoupling both. 

Second, we want to stop the collapse along any of the prin- 
cipal axes according to a physically motivated virialisation con- 
dition. Virialisation must be invoked to prevent the axes from 
collapsing to zero, and thus to be able to follow the entire 
collapse of an ellipsoid, i.e. the collapse of its three principal 
axes. Conventionally, the collapse of each axis is stopped when 
a,(a) - Q.nia, where a, is the scale factor of the i-th axis and a 
is the background scale factor dBond & MversI 19*961 ; I Sheth et al.l 
2001). In this way, A - 0^/(010203) - 178 for spherical collapse 
in an Einstein-de Sitter universe at the time when the third axis 
virialises. However, the value 0.177 has no fundamental physi- 
cal motivation, and there is no guarantee for it not to be different 
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for ellipsoidal rather than spherical collapse, or when cosmolo- 
gies other than EdS are to be considered. We introduce a general 
virialisation condition based on the tensor virial theorem that 
avoids introducting such an uncalibrated factor. We find substan- 
tial changes on both Ay and 6^ from both modifications and point 
out several discrepancies of our results with earlier studies. 

We introduce the ellipsoidal-collapse model including our 
modifications in Sect. |2l present our results in Sect. [3] and con- 
clude in Sect.m 

2. The model 

In this Section, we shall briefly review the ellipsoidal-collapse 
model of Bond & Mveri (Il996h for cosmologies with a cosmo- 
logical constant introducing the scale factor a as time variable. 
Furthermore, we shall present a physically motivated virialisa- 
tion condition to stop the collapse of each axis, and show how to 
find the proper initial ellipticity and prolaticity as a function of 
mass and virialisation redshift. 



2.1. The evolution equations 

Let fl, - Ri/Rpk be the dimension-less principal axes of the el- 
lipsoid, where Ri with 1 < / < 3 are its dimensional semi-major 
axes, and Rpi^ the size of a spherical top-hat corresponding to a 
mass M = (4n/3)phR^^ with the cosmological background den- 
sity pb- The evolution of the three principal axes fl, with time f 
in a cosmology with a cosmological constant A = (87tG/c^)pA 
is given by 



d^fl, 8 (1 6 bi 

= -nGp^Oi - 47TGpbfl; I - + 3 + + '^^'^t.'' 



(2) 



dBond&MversI 119961) . where G is the gravitational constant 
and c the speed of light. The density contrast of the ellipsoid 
with respect to the background density v& 6 - (p - Pb)lpb - 
fl-'/(flifl2fl3) - 1. The parameters /?, and A^^ij denote the internal 
and external contributions to the gravitational tidal shear which 
occur because of the deviation from sphericity. Generally, the 
total tidal field is described by the tidal field tensor T with the 
elements Tjj - 5^<l>p/(5x,c)xj) - Tintjj + Tf^-^^jj, where Op de- 
notes the peculiar gravitational potential, and Tint,,/ and Text,,) are 
the internal and external contributions to the shear, respectively. 
After a transformation into the ellipsoid's eigensystem, which is 
the same as the eigensystem of T in this model, the internal shear 
can be evaluated as 



dr 

bi(t) = ai(t)a2(t)a3(t) — — — — 
Jo [fl2(f) + l]nLi[«'(0 + 

while the external shear can be approximated by 



l]i/2 3 



Aext,i(t) = 



Mto) - 



(linear approx.), 
(non-linear approx.) , 



(3) 



(4) 



where D+ is the linear growth factor, and the /I, are the eigenval- 
ue softheZeFdoyich deformation tensor See the first appendix 
of lBond & MversI ( 1996) for details of the calculation. Here and 
in the following the index '0' refers to initial values. 

In linear approximation, the environment into which the el- 
lipsoid is embedded evolves completely independently of it, 
whereas in the non-linear approximation it is tightly coupled to 



the ellipsoid. In Sect. l3.1l we shall introduce the hybrid model as 
a third approximation for the external shear. 

We can rewrite Eq. (|2]i by using the scale factor a as time 
variable using Friedmann's equation cp- - HlE\a), where Hq 
is Hubble's constant. The expansion function of the universe, 
^(fl), introduces the dimension-less density parameters of matter 
and dark energy today, Q.^ and Qa, respectively. This gives 



d^ 



1 E'(a) 
a E(a) 



dfl, 
dfl 



30, 



2fl5£2(fl) 



C,(fl) ■ 



Qa 



fl2£2(a) 



(5) 

where a prime denotes differentiation with respect to a, and C,- = 
(1 +5)/3 +/7,/2 + /lexu- Equation (|5j defines a set of three coupled 
second-order differential equations, for which we need six in- 
dependent initial conditions compatible with the Zel'dovich ap- 
proximation for early times. These are provided by 



fl,(flo) 

dfl, 
dfl 



:fl()[l - /i,(flo)] . 



= 1 - /l;(flo) - 



dlnD. 



dlnfl 



/l,(flo) 



(6) 

l-2^,(flo), (7) 



since D+{a) « a for fl <K 1. Choosing flo = 2 x 10"'', this is 
comfortably fulfilled. In Ap pendix lAl we co mpare Eqs. (|5}Q 
with the results presented by iMonacol (Il997h and see that they 
diff^er. 

In the following, we shall assume that the eigenvalues A, of 
the Zel'dovich tensor are ordered as Ai > A2> A^, which implies 
that the ellipsoid first collapses along the direction 1, fl] — > 
first. At that time (5 — > oo, and the collapse of the remaining 
two axes can no longer be followed so that we have to add a 
virialisation condition for each axis preventing their collapse to 



2.2. Tine viriaiisation condition 

Conventionally, the collapse of each axis is stopped when fl,(a) = 
0.177fl, where fl,- is the scale factor of the i-th axis an d a the 
background scale factor (e.g. Bond & Mvers 1996; Shet h et aP 
I200lh . In this way, A - fl-'/(flifl2fl3) = 178 for spherical col- 
lapse in an Einstein-de Sitter (EdS) universe at the time when 
the third axis is assumed to virialise. However, the value 0.177 
has no fundamental motivation in the physics of the collapse, 
and there is no guarantee that it remain unchanged in the case of 
ellipsoidal instead of spherical collapse, or in other cosmologies 
than Einstein-de Sitter 

We shall thus follow a dififerent approach and present a phys- 
ically well-motivated virialisation condition in the following to 
stop the collapse of ea ch axis individually. We sta rt from the ten- 
sor virial theorem (see lBinnev & Tremaind 19871 p. 213, p. 280), 



1 d^Ii, 



2 df2 



2Kij + Hi J + Wij + Vij 



(8) 



where I is the moment of inertia tensor, K and fl are the con- 
tributions to the kinetic energy tensor coming from ordered and 
random motions, respectively, W is the potential energy tensor, 
and V is the external potential energy tensor. Their elements are 
generally defined as 



Iij= dxpxiXj, Wij = - dxpxi—, 
J'V Jv o^j 

Kij =\ ^ d^xpvivj , n,j = J d^xpojj , 



(9) 
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and 



1 / 5cDext 

= dxp\xi— — 

2 J<v \ dxj 



+ Xi 



dxi 



(10) 



where p is the density of the fluid, <b and Oext are the gravi- 
tational potentials of the ellipsoid itself and its surroundings, 
respectively, 'V is the volume which is integrated over, cr]. = 
ViVj — vivj are the velocity dispersions, and the bar indicates av- 
eraging over 'V. 

We now specialise to the case of a homogeneous ellipsoid. 
For a stable mass configuration, the left-hand side of Eq. ([8]) has 
to vanish for each component of the inertial tensor individually. 
Since the ellipsoid is assumed to be at rest and the ellipsoid's 
eigensystem is chosen as a reference frame, v, = so that Kij = 
and Tlij - J^d^xp{vj)6jj, with the Kronecker symbol 6jj. Note 
again that in this framework the eigensystems of the overdense 
ellipsoid and the gravitational tidal field are identical. For a ho- 
mogeneous ellipsoid v',(x,) = (fl,7fl,)x,-, thus 



(11) 



w. 



with the mass M of the ellipsoid. The sum Wtj + Vij can be eval- 
uated using -V(<1> -I- Oext) - X and Eq. (O to be 

r , d^UiXi 1 , /87tG \ 

(12) 

Requiring that Eq. dHJ is fulfilled for each axis separately and 
introducing the scale factor a as time variable yields the viriali- 
sation conditions for the three axes a,-. 



1 



3Q„ 



C, - Q/ 



(13) 



When this condition is fulfilled for an axis together with a, < 0, 
its collapse is stopped and its size is frozen in. 

We emphasise that the former equation is consistent with 
the virialisation condition for spherical collapse in the EdS uni- 
verse, where Ry/Rta = 0.5. The subscripts 'v' and 'ta' denote 
the time of virialisation and turn-around, respectively. For the 
special case of EdS, Eq. ( fT3] ) becomes 



R 



GM 



with a\ - a2 - ut, = R. Using the parametric solution 



R 



Rq^ 

2do 



(1 -COS0), 



3fo 



(6 - sin 0) 



(14) 



(15) 



([Engineer et al.ll2000l) together with the relations Rq/So - Rm 
and Rl = 9GMtl/2 indeed gives the expected result that Eq. (fT4l l 
is satisfied when R - /?ta/2. 

For bound objects in the ACDM model, the virialisation con- 
dition is 

<£kin) = -^<£pot) + <£A>, (16) 

where (/Skin) is the average kinetic energy of the halo, (fipot) 
its average potential energy, and {E^} the average effective po- 
tential energy contributed by the cosmological constant. Using 
(Epot) = -3GM^/(5R) and (^a) = -AMR^/IQ for the homo- 
geneous sphere as well as energy conservation between turn- 
around and virialisation, one arrives at a cubic equation for 

Ry/Rm, 



(2 + /,)|^ + l 



0, 



(17) 



at 




tensor virial theorem, z=Z| 
tensor virial tineorem, z=z. 



normal virial tiieorem 
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Fig. 1. Comparison of the virialisation conditions from the scalar 
and tensor virial theorems. The curves differ slightly because 
they refer to different times. While the tensor virial theorem was 
evaluated at virialisation, the scaler virial theorem was evaluated 
at collapse. If we choose the collapse of the dark-matter halo as 
reference for z, both yield exactly the same result. 



with T] = ARl/(3GM) dLahav et al.ll 19911) . The relevant solution 
of this equation agrees precisely with the condition derived from 
the tensor virial theorem, Eq. ( flj] ), as Fig.[T]shows. The small de- 
viation occurs because we choose the time of virialisation rather 
than collapse as a reference (z = Zy) when using the tensor virial 
theorem. If we use the collapse time (z = Zcoi) instead, both con- 
ditions yield identical results. 

2.3. The initial ellipticity and prolaticity 

Equations (|6] [T]) imply that one has to choose initial values for 
the Zel'dovich deformation tensor to define the initial deviation 
of the principal axes and their time derivatives from the back- 
ground. We shall show how they are chosen appropriately so that 
they comply with the assumed Gaussian nature of the Universe's 
initial conditions and represent a statistical average of haloes 
with the same mass M and virialisation redshift Zv but different 
shapes. 

Starting from the probability distribution for the eigenval- 
ues of the Zel'dovich deformation tensor, Eq. ([1]), the condi- 
tional probability distribution for the ellipticity e > and the 
prolaticity -e < p < e, defined as e = (Ai - Ai,)/( 26) and 
p_=_(Ai-2A2+Ai)/(26), respectively, was derived bv iSheth et al.1 
hoOlh to be 



g(e,p\S) 



1125 



(18) 

To compute a statistical average of any quantity ^{e,p) for a 
given halo that depends on e and p, one should marginalise over 
the distribution g(e, p), 



Jo J-e 



dp^ie,p)g(e,p) . 



(19) 



However, a Taylor expansion of (^) up to second order in e and 
p around their mean values (e) and {p} under the distribution 
g(e,p) gives 



1 5V 

(^)=^((e>,(p)) + - ^ 



<p>,(p> 



2^1 
"^^2 5? 



(20) 



(e),<p) 
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where cr^ and cr^ are the variances for e and p according to the 
distribution g. Up to first order, (^) = ^{{e),{p)). Deviations 
occur only at second order The expectation values for e and p 
as well as their variances are given by 



{e) 



3cr 



(197t-54)o-2 
607t(52 



2 ^ 



(21) 
(22) 



Since the variances are oc (cr/6)^, and o-/6 < 1, using the approx- 
imation {^) X! f ((e), {p}) introduces only a small error, which is 
-1% for 4 and -3% for Av for ACDM and OCDM. For EdS, 
the error is larger and ~3% for 6c and ~10% for Ay. But since 
the latter cosmology is scientifically only of low relevance and 
usually serves as a reference model only, the usage of the for- 
mer approximation is well justified. Instead of sampling ^(e, p) 
at several points for e and p, one only has to evaluate it once for 
(e) and {p} = 0. 

Generally, the eigenvalues /I, are related to e and p by (see 
e.g. iBond & Mvers..l996: Bardeen et al...l986,) 



6 6 a 

Ai^ -(\+3e+p)^ - + —= , 
3 3 VlOTt 

6 6 
A.^-(l-2p)^-, 

6 60- 
3 3 VlOTt 

where we have set e - {e) and /? = in the last step. 



(23) 
(24) 
(25) 



3. Results 

In this Section, we show the results of the ellipsoidal-collapse 
model for the parameters 6^ and Ay for three different cosmolo- 
gies and discuss how they are affected by the choice of the 
external-shear model. 



3. 1. Influence of the external shear 

Figure |2] shows the evolution of the three principal axes and the 
eigenvalues of the external shear for a reference ACDM model 
with - 0.3, Qa = 0.7, and erg = 0.8 for three different 
models of the external shear and e - (e), p - Q. The shear in 
the linear approximation evolves completely smoothly over the 
entire collapse time. The evolution of each eigenvalue is given 
by D+(a) so that according to Eq. /lext.i > 0, /lext,2 = 0, and 
Aext.3 - -'iext.i < at all timcs. 

This is different in the non-linear approximation: At early 
times, the evolution of the A^xij is the same, but soon thereafter 
they start evolving non-linearly and /lext,2 becomes slightly neg- 
ative. Noticeably there are steps in the evolution of the exter- 
nal shear whenever an axis virialises because A^xtj flifl2fl3, 
and the evolution of this volume factor changes after virialisa- 
tion of each axis. In the right panel, we introduce the hybrid 
approximation: Initially, the evolution of /text,/ is described by 
the non-linear model. When one of the axes turns around, how- 
ever, the corresponding eigenvalue of the external shear contin- 
ues evolving linearly, i.e. its value at turn-around is then scaled 
by D+(fl)/D+(fl,a). 

We believe that the hybrid model is the preferred model for 
the evolution of the external shear since it takes into account 
that the evolution of the ellipsoid itself and its vicinity should be 
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Fig. 3. 6c (left panels) and Ay (right panels) for three different 
models of the external shear Both quantities were calculated for 
a 10'^ Ma/h dark-matter halo in the reference ACDM model. 
Top panels: e = (e) as a function of virialisation redshift Zv 
Bottom panels: e - 0.2 for all Zv In both cases, p - 0. 



tightly coupled in the beginning. At turn-around, however, they 
are definitely decoupled so that choosing this moment to switch 
from non-linear to linear evolution seems appropriate. Hence, 
we will use the hybrid model for the evolution of the external 
shear in the following. 

The influence of the external-shear model on the parameters 
6c and Ay is shown in Fig.|3] For a given mass and a given virial- 
isation redshift Zy, the initial overdensity is chosen such that the 
third axis of the ellipsoid finally virialises at Zy. Both parameters 
can then be calculated by 



D+(av) 
D+iao) 



J] Mao) , 
1=1 



nm(flv)fly 



fll(flv)a2(flv)fl3(flv) 



(26) 



In this case. Ay is defined with respect to the critical density as 
it can be primarily found in the literature. For the definition with 
respect to the background density one simply has to omit Q^iay). 

Using (e) as a function of virialisation redshift and p - (top 
panels of Fig.O, the dependence of 6c and Ay on Zy is almost the 
same for the non-linear and the hybrid models. However, both 
differ from the linear model, showing that the external shear 
is most important at the beginning of the ellipsoid's evolution. 
At that time, the non-linear and the hybrid models agree. While 
6c(Zv) is always smaller in the linear compared to the other two 
models, the curves for Ay(zy) cross. This reflects the circum- 
stance that the initial overdensity in the linear model is different 
from that in the two other models, leading to a different initial 
ellipticity and therefore to a completely different evolution his- 
tory. This can be seen in the bottom panels of Fig. [3] for which 
we have chosen e = 0.2 independently of Zy. They also clearly 
show that a varying initial ellipticity drives primarily the evo- 
lution of 6c, whereas Ay also strongly varies for fixed e. In this 
case, both 6c and Ay are smaller in the linear-shear model com- 
pared to the non-linear and the hybrid models. 

3.2. Parameters as function of mass and redshift 

Before we present general results for the parameters 6c and 
Ay, we should comment on a subtle but very important issue: 
Whenever we want to compare our results with the ordinary 
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Fig. 2. Evolution of the principal axes and the external shear for a 10 M0//1 dark-matter halo with e - (e) and p = in the reference 
ACDM model. Left panel: Linear approximation. Central panel: Non-linear approximation. Right panel: Hybrid approximation. 
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Fig. 4. 6z (left panels) and Ay (right panels) for three differ- 
ent cosmological models and fixed initial e and p. Top panels: 
Spherical case (e - p - 0). Bottom panels: Triaxial case with 
e = 0.2 and p = -0.1. 

spherical-collapse model, we have to keep in mind that we cal- 
culate all quantities at the time when the third axis virialises. 
Thus, we also have to compare these quantities those from the 
spherical-collapse model that are also calculated at the time of 
virialisation and not of collapse, i.e. when R - /?ta/2 and not 
R = for EdS. This leads to slightly lower reference values of 6^ 
and Ay since Zcoi < Zv Here the s ubscript 'col' den otes collapse. 
Using the parametric solutions of lOhta et al.l (|2004|) for the linear 
and the non-linear overdensity, 

3 [3 f'^ 9 (6* -sin 61)2 

= - -(6 -sine) , A„i = -i (27) 

5 [4 J 2 (1 - cos 6')-' 

respectively, and 6 = 37T/2 at viriaUsation, we find that 6c - 
1.583 and Ay = 147 for the EdS universe at R - R^ indepen- 
dent of Zy. Recently, Lee & Ng , C2009.) arrived at the same values 
when accounting for the time of virialisation instead of collapse. 

The top panels of Fig. |4] show 6^ and Ay for three different 
cosmologies for e - p - Q, i.e. spherical systems. The OCDM 
cosmology is the same as our reference ACDM model except 
that Oa = 0. For the EdS model we set = 1 and Oa = 0. 
Indeed, for the EdS universe the constant values derived analyt- 
ically are also reproduced by solving Eq. (|5]l numerically. This 
demonstrates again that Eqs. dSlflJll are fully consistent with the 
well-known spherical-collapse model. Note that the qualitative 
behaviour of Av(zv) is the same as Av(Z(;oi) (compare e.g. with 



Fig. 5. 5c (left panels) and Ay (right panels) as a function of halo 
mass and virialisation redshift for three different cosmologies 
marginalised over e and p. Top panels: Dependence on virialisa- 
tion redshift Zy for a halo with M - lO''* Mq/Zi. Bottom panels: 
Dependence on mass M for a halo with Zy = 0. Thin dashed 
lines show the reference values from spherical collapse in the 
EdS universe, 6^ = L583 and Ay = 147. 

iBartelmann et al.ll2006l) . However, there is a difference for the 
critical linear overdensity whose shape as a function of Zy differs 
substantially from the shape as a function of Zcoi- This should il- 
lustrate that the time chosen in the model when virialisation ac- 
tually occurs (zy or Zcoi) can already have substantial impact on 
the qualitative behaviour of relevant quantities as a function of 
redshift, hence it is not necessarily a consequence of ellipsoidal 
collapse alone. 

The bottom panels of Fig.|4]show 6^ and Ay for a triaxial halo 
with e - 0.2 and p - -0. 1 . Also in this case both parameters 
are independent of Zy for the EdS universe, although 6^ changes 
from 1.583 to 2.058, while Ay stays almost the same, 148 instead 
of 147. Interestingly, the most drastic changes in the shapes of 
both parameters occur for the OCDM model, for which the total 
density is only approximately a third of the critical density, while 
for the ACDM model their changes are small. 

The influence of the halo mass M and the virialisation red- 
shift Zy on both (Jc and Ay are illustrated in Fig. |5] For both 
decreasing mass and decreasing redshift, 6^ is a monotonically 
decreasing function, approaching the reference values from the 
spherical-collapse model in the EdS universe for large masses 
and high virialisation redshifts, the situation for Ay is much more 
complicated. For all three models, it has a minimum at redshifts 
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Fig. 6. Comparison of the fitting formula of ISheth et al.l (1200 lb 
for 6c{M, z) with the results of this work for the reference ACDM 
model. Left panel: 6c as a function of mass and z = 0. Right 
panel: 5^ as a function of redshift and M - 10^"* Mq/Zi. 

4-5 and at a mass of ~ lO'^ Mq/H. For smaller values of mass 
and redshift, it is a monotonically decreasing function of both M 
and Zv for the OCDM and EdS models. It is also monotonically 
decreasing as a function of M in the ACDM model, but reaches a 
maximum at Zv ~ 1 ■ This is a direct result of the definition of Ay 
with respect to the critical density. If it was defined with respect 
to the background density, the factor Qn,(flv) would not appear 
in Eq. ( |26] |, and all three curves would increase with decreasing 

Zv. 

In Fig.|5]one can clearly see that the intervals that are covered 
for both 6c and Av are the largest for the EdS model, indicating 
that there is a stronger dependence of the ellipsoid's evolution 
on the total amount of matter in the Universe compared to the 
size of the cosmological constant if varying initial elUpticities 
are taken into account. 

For either M — > oo or Zv — > °o, both 6c and Ay must reach 
the reference values for spherical collapse in the EdS universe 
since the initial ellipticity (e) cc cr/6 and cr both decrease with 
increasing mass, and 6 has to be higher the earlier the structure 
is required to collapse. This expected behaviour can be clearly 
seen for 6c, whereas for Ay this happens finally for very large 
M and Zy. We should stress again in this context that a crucial 
portion of the dependence on mass and virialisation redshift is 
driven by the change in the initial ellipticity, comparing Figs.|4] 
and ID 

We compare in Fig.|6]the results of this work for 6c as a func- 
tion of mass and redshift with the fitting formula of ISheth et al.l 
(1200 lb given by 



6ciM,z)^6c,spb(z)\l+0A7 



o-2(M,z) 



<.sph(^) 



{).615^ 



(28) 



where i5c,sph(z) is the redshift-dependent linear overdensity of the 
usual spherical-collapse model. The left panel shows that the de- 
pendence on mass is similar for both the fittin g formula and the 
result of this work. However, the fit bv .Sheth et all (l2001b is less 
steep as a function of mass resulting from the diff'erences of their 
underlying work to ours: First, they use the linear-shear model 
instead of the hybrid model, leading to a larger 6c for all redshifts 
as shown in Fig. [3] Second, the collapse of each axis is stopped 
using the artificial condition a, = Q.llla, leading to 6c - 1.686 
as a reference value at high mass. Thus, the fitting formula pro- 
vides a larger value for large masses, while our approach leads 
to a (5c — 1.583 for M — > oo due to the virialisatio n condition 
that we apply. Third, for a given initial overdensity ^. ISheth et alj 
(12001") use the most probable value emp = (c/6)/ Vs instead of 
the expectation value (e) given by Eq. (l2Tl i. which leads to initial 
ellipticities that are slightly too low so that the asymptotic Umit 
for high redshifts is reached earlier than in our case. 
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Fig. 7. Influence of the initial ellipticity e and prolaticity p on the 
parameters 6c {top panels) and Ay (bottom panels) for three dif- 
ferent cosmological models. Left panels: ACDM. Central pan- 
els: OCDM. Right panels: EdS. The large circle in each figure 
indicates the value of the respective parameter for e - (e) and 
p = 0. The small circles indicate values for a combination of 
e € [<e> + (Te/2, <e) + crj and p e [±e/4, +e/2] (see Eqs.lITI 
I22I 1. For a fixed e both 6c and Ay grow for decreasing p. 



The right panel of Fig. |6] shows a similar behaviour for 6c 
as a function of redshift. The differences between the fitting for- 
mula and the result of our work again occur due to the different 
virialisation time and condition, and the difference between the 
most probable and the expectation value of e. 

For applications, fitting formulae for both 6c and Ay may 
be useful. We provide expre ssions here which are inspired by 
Eq. dUl) of lShefli et all (l200ll) . For 6ciM, z) we suggest 



(Jc(M,z) = 5c,sph(Zv)P +^ 



0-2(M,Zy) 



^c.sph(^v) 



(29) 



where b - 0.6536, c = 0.6387, and both <5c,sph and cr^ are 
cosmology-dependent quantities. Note again that the redshift of 
virialisation, Zv, has to be chosen as reference for z for both ^c.sph 
and cr^. For a spatially-flat ACDM model with matter density pa- 
rameter in the range e [0.2,0.4], M e [10", 10'^] Mq h-\ 
and z e [0, 10], the maximal error is ~1.8% with a mean error of 
-0.4%. 

A similar functional dependence can be found for Ay(M, z). 
Only a small correction term has to be added to arrive at a satis- 
factory accuracy. We find 

Ay(M,z) = Ay.,ph(Zy) [a + bcr^'(M,Zy) + d{l+ z^f'^ \og'^\M)\ 

(30) 

with a = 0.3819, b = 0.5379, c = 0.7589, and d = 3.456 x lO^^. 
In the same range as above, the maximal error is ~5% with a 
mean error of 

3.3. Influence of initial ellipticity and prolaticity 

In Fig.|2]we plot both 6c and Ay as a function of the initial ellip- 
ticity e and prolaticity p centered around their expectation values 
given by Eq. (1211 1 for three different cosmologies. For increasing 
e and decreasing p, both paramet ers grow qualitativ ely in the 
same way as already reported by ISheth et al.l (1200 ll) (cf. their 
Fig. 1). Quantitative deviations arise from the differences in the 
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applied algorithm as discussed in Sect. 13. 21 For a given mass and 
virialisation redshift, the initial overdensity for the EdS universe 
is larger compared to both ACDM and OCDM due to a shorter 
physical time interval that corresponds to the same redshift inter- 
val, resulting in a larger (e) and cr^, but also in larger curvatures 
of 6c and Ay with respect to e and p. These are the sources of 
the larger error in the approximation (^) ^ ^((e), {p)) discussed 
in Sect. 12.31 Since the redshift-time relation is not very different 
between ACDM and OCDM, the dependences of 6c and Av on e 
and p are comparable. 

3.4. Mass function 

Using Eq. ( |29] ). we are able to construct the mass function of 
dark-matter haloes using the extended Press-Schechter formal- 
ism developed by Bond et al. ( 1991) and Lacey & Cole ( 1993]), 
which is based on the first-upcrossing distribution of the density 
contrast 5 as a function of the "time variable" S = cr^ jM). We 
shall proceed similarly as lSheth & Tormeril ( ll999ll2002h and de- 
fine the scaled variable v = 6^ i./S to derive the mass function 

c.sph' 

for ou r standard ACDM cosmo logy. 

As lSheth & TormenI (l2002h pointed out, expressing the first- 
upcrossing distribution / as a function of v has the advantage 
that it is only necessary to calculate /(v) for a barrier of height 
B{v, z) at one arbitrary redshift to infer the mass function n(M) 
at any other redshift by a simple rescaling. For a given first- 
upcrossing distribution /(v), the differential mass function can 
be calculated using the relation 



^ ^ M^dlnM ■'^ ^' 



(31) 



where pb is the background density of the Universe. 

First, we want to find an accurate fit to the first-upcrossing 
distribution of a moving barrier which is given by the mass- 
dependent linear overdensity parameter of the ellipsoidal col- 
lapse, 

B(v) = (5c,sph (l + 0.6536 y-°-^^**') (32) 

(see Eq. |29] ). The parameter i5c,sph is evaluated at Zv - 0. We 
ran one million random walks and recorded the first-upcrossing 
values for v e [0.01,20] in 100 equidistant bins in logarithmic 
space. The resulting distribution v/(v) is nicely expressed by the 
function 



v/(y) = A [1 -H iav)-P] 



■ exp 



av B(v) 



2 6, 



'c,sph 



(33) 



Thus, our suggeste d fitting formula is a mix tur e of the functional 
forms proposed bv lSheth & TormerJ d 19991) and lSheth & TormenI 
(|2002|) . The remaining best-fit parameters are A - 0.357, p - 
0.212 and a - 1.171. The result is shown in Fig.|8] 

Second, to find a viable mass f unction from th e first- 
upc rossing distr i bution , we proceed as ISheth & TormenI d 19991) 
and ISheth et al] (1200 lb . normalise /(v) to unity and rescale the 
variable a such that we are in agreement with the standard 
[Sheth & Tormen mass function and a mass function based on 
A^-body simulations proposed bv ICourtin et aki (120 lOl) . The lat- 
ter is based on a first-upcrossing d istribution that has the s ame 
functional form as that proposed bv lSheth & TormenI (1 19991) . but 
slightly different best-fit parameters. 



y/coui1in 



(v) = A [l + (av)-'''] y^exp(-y) , (34) 



0.01 




0.001 



0.0001 - 



'csph'' 



Fig. 8. Comparison between the first-upcrossing distribution of 
the moving barrier, Eq. ( l32l i. inferred from an ensemble of one 
million random walks, and the fitting formula, Eq. ( 1331 ). 



Sheth & Tormen (2001) 
Courtinet al. (2010) 
this work 




Fig. 9. Comparison of the I Sheth & TormenI and the lCourtin et al.l 
mass functions with the mass function derived from the rescaled 
upcrossing distribution, Eq. ( |33] ). based on our treatment of the 
elUpsoidal-coUapse dynamics. 



with A = 0.348, a = 0.695, and p = 0.1. Note that in their 
definition of v, the linear density contrast (Jc.sph has to be taken 
at collapse. Normalising the first-upcrossing distribution based 
on the moving barrier of our ellipsoidal-collapse model to unity 
yields a rescaled parameter A ^ A' = 1 .364 A. We compare the 
resulting mass fun ction with those bv ISheth & Torme nl ([T99l 
and Courtin et al.l QOIO) for three different redshifts in Fig. |9] 
The par ameter a was resc aled by a — > a' = 0.625 a. Deviations 
from the lSheth & Torme r?mass function at high mas ses oc cur at 
large redshifts which is compatible with the Courtin 'et alJ mass 
function. Overall, our proposed mass function lies in between 
these two, s uggesting that diff erences between A^-body simula- 
tions and the lSheth & TormenI mass function might be due to an 
imprecise treatment of the ellipsoidal-collapse dynamics. 



4. Conclusions 

We have reconsidered the collapse of a homogeneous, triaxial 
ellipsoid in an expanding backgro und universe and extended the 
treatment of iBond & MversI (Il996<) in two ways: 
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- We have introduced a physically motivated criterion for the 
onset of virialisation along each principal axis of a collaps- 
ing ellipsoid. We derive this criterion from the tensor virial 
theorem, demanding that the inertial tensor stabilises and its 
second time derivative vanishes. This approach is a general- 
isation of the usual virialisation condition for the spherical- 
collapse scenario and thus fully consistent with it. It replaces 
the conventional requirement that virialisation is assumed 
along an axis / when the dimension-less semi-major axis a,- 
of the ellipsoid along that axis reaches 0,(0) = QAlla. 

- We have introduced a hybrid model for the influence of 
the external gravitational shear field acting on the ellipsoid, 
supple menting the two diff'erent models of iBond & MversI 
(Il996l) . The principal axes of the shear field evolve with 
those of the halo before it turns around and decouples from 
the background expansion, and then continues to evolve lin- 
early while the halo evolves non-linearly and collapses. We 
have shown that the differences between the hybrid and the 
non-Unear model are relatively small. 

For a given initial ellipticity e and prolaticity p, and for a 
specified virialisation redshift Zv, the ellipsoidal-collapse model 
then gives a unique answer for the linear density contrast 6c at 
virialisation, as well as for the overdensity Ay at that time. The 
probability distribution for e and p, conditional on the density 
contrast 6, is determined by the probability distribution of the 
eigenvalues of th e Zel'dovich deformation tensor, as shown by 
[She57 et al.l ((2001). It is characterised by variance of the matter- 
density fluctuations on a scale fixed by the halo mass required. 
We have shown that the marginalisation over e and p can be 
simplified by evaluating the ellipsoidal collapse at their mean 
values (e) and {p} = 0. 

Our main results are as follows: 



- The collapse parameters 6^ and Ay depend only weakly on 
the model for the external gravitatio nal shear. The hybrid 
model and the non-linear model by iBond & MversI (Il996h 
give approximately the same results. 

- When supplied with our virialisation condition derived from 
the tensor virial theorem, the ellipsoidal-collapse model re- 
turns values (5c and Ay that differ substantially from those ob- 
tained with the spherical-collapse model. Depending on halo 
mass and redshift, deviations of order (20 . . .50)% are com- 
mon. After marginalisation over e and p, 6c and Av increase 
with decreasing halo mass and with decreasing virialisation 
redshift. 

- Both parameters increase with increasing initial ellipticity e 
and decr ease with increasin g prolaticity p, as already sug- 
gested bv lSheth et all (l200lh . 

- Our results for 6c qualitatively confirm the dependence on 
ha lo mass and virialis ation redshift given by a fitting formula 
bv lSheflietan(l2001h . Deviations in particular at low mass 
and high redshift occur due to the differences in the virial- 
isation condition, the model for the external shear, and the 
marginalisation over e. 

- The mass function based on our refined treatment of the 
ellipsoidal-c ollapse model is in good agreement with those 



proposed by 



ICourtin et al. 



Sheth & Torineiil (Il999h for low redshifts and 
(120 10^ fo r high redshifts. This suggests that 



difl^erences between the'Shet h & TormenI mass function and 
results from A^-body simulations at large redshifts may oc- 
cur due to an imprecise treatment of the ellipsoidal-collapse 
dynamics. 



We now plan to introd uce the ellipsoidal-collapse mod el into 
the approach presented bv 'Angrick & BartelmannI ( l2009l) to im- 
prove the prediction of the X-ray temperature function for galaxy 
clusters based on the gravitational potential rather than on the 
density contrast. 
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Appendix A: Comparison to a previous study 

Here we compare our results for the evoluti on equat i ons of the 
collapsing ellipsoid with those presented bv lMonac ol ([T997l) for 
a flat ACDM and an OCDM model. 

Starting from Eq. (|5]l, we can replace E(a) and E'(a) for a 
flat ACDM model using 

E(a) = V"mfl-3 + (1 - S^m) , = — (A.l) 

2E(a) 

since Qa = 1 - Qm and the curvature parameter Qk = 0. This 

gives 

d^Gj _ l-2(Q„'-l)a^ da,- ^ 3Ci-2a\£l^^ - 1) ^ _ ^ 
dfl^ 2fl [1 + (0„/ - l)fl3] dfl 2fl2 [1 + (Q^i - l)fl3] 

(A.2) 

This equation diff'ers from Eq. (Bll) of iMonacol (Il997h in the 
second and third term: a factor a- in the denominator of both 
terms was omitted. Additionally, the vacuum term oc as well 
as a factor 3 were not included in the nominator of the third term. 
For an OCDM model, we have 



E(a) = V^^mfl"^ +(1 -Om)fl"^ 



E'(a) : 



2E{a) 



(A3) 
(A.4) 



since Qa = and Qk = (1 - f2m)fl ^. Inserting this again into 
Eq. (|5]l yields 

d^ai ( r „ 1 ni-i da,- 
^-{24l.(0--l)«]}^^ 

+3 {2fl2 [1 + (0„' - l)a]] ' C/fl/ = . 



Equation (B12) of lMonacol (1 1 997h is again slightly diff'erent: The 
factor 3 in the last term was omitted. 

There is one last diff erence concernin g the initial conditions: 
Comparing Eq. (B17) of lMonacol (1 1997h with Eq. (jgl of this pa- 
per, one can find an additional factor qq in front of /l,(ao) which 
should be dropped. 
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